Ginzburg-Landau theory for the conical cycloid state in multiferroics: applications to CoCr204 
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We show that the cycloidal magnetic order of a multiferroic can arise in the absence of spin and lattice 
anisotropics, for e.g., in a cubic material, and this explains the occurrence of such a state in CoCr204. We 
discuss the case when this order coexists with ferromagnetism in a so called 'conical cycloid' state, and show 
that a direct transition to this state from the ferromagnet is necessarily first order. On quite general grounds, 
the reversal of the direction of the uniform magnetization in this state can lead to the reversal of the electric 
polarization as well, without the need to invoke 'toroidal moment' as the order parameter. 
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I. INTRODUCTION 

Ferromagnetism and ferroelectricity are two of the most 
well-known and technologically relevant types of long range 
ordering that can occur in solids. It is therefore of paramount 
interest and importance that in a class of ternary oxides, 
known as "multiferroics", both types of order seem to co- 
exist with the possibility of interplay between long range 
magnetism and long range electric polarization'-^'^-^. The 
recently discovered new class of multiferroics with strong 
magnetoelectric effects often display the coexistence of a 
spatially modulated magnetic order, called 'cycloidal' or- 
der, and uniform polarization (P), which is induced by the 
broken inversion symmetry due to the modulation of the 
magnetization^'^. Since P is inherently of magnetic origin, 
unusual magnetoelectric effects, as displayed by the ability 
to tune the polarization by a magnetic field which acts on 
the cycloidal order parameter, are possible, opening up many 
application a^i^i^'^i ' i- ' 2, i3, 14, 1 5 ^ Among this exciting class of 
materials, the cubic spinel oxide CoCr204 is even more un- 
usual, since it displays not only a non-zero P and a spatially 
modulated magnetic order, but also a uniform magnetization'^ 
(M) in a so-called 'conical cycloid' state (see below). The 
uniform component of M provides an extra handle^ with 
which to tune P, as has been recently demonstrated'^. The 
low value of the required tuning magnetic field ^ .5 T, makes 
this material even more experimentally appealing. 

The ability to tune P by tuning the uniform part of M poses 
a theoretical puzzle, since, in existing theories, the uniform 
piece of M should not influence the polarization at all^j^J^ill. 
This has lead to the introduction of the 'toroidal moment', 
T = P X M, as the real order parameter characterizing the 
conical cycloid state of CoCr204'^. In this Letter, we explain 
this unique phenomenon and the other interesting aspects of 
the physics of the conical cycloid state by developing a phe- 
nomenological Ginzburg-Landau (GL) theory. Additionally, 
the rotationally invariant form of the theory proves that both 
the ordinary and the conical cycloidal orders, with the result- 
ing multiferroicity, are possible even in systems without easy 
plane spin and easy axis lattice anisotropics. This is impor- 
tant since earlier models^sliiS of the cycloidal state depend cru- 



cially on such anisotropies. However, such anisotropic mod- 
els can not explain the presence of the cycloidal state in cubic 
systems like CoCr204, where such phases are also observed 
despite the fact that their cubic symmetry forbids such easy 
plane and easy axis anisotropies. 

CoCr204, with the lattice structure of a cubic spinel, en- 
ters into a state with a uniform magnetization at a temperature 
Tjn — 93 K. Microscopically, the magnetization is of ferri- 
magnetic origin'^, and in what follows we will only consider 
the ferromagnetic component, M, of the magnetization of a 
ferrimagnet. At a lower critical temperature, Tc = 26 K, the 
system develops a special helical modulation of the magneti- 
zation in a plane transverse to the large uniform component. 
Such a state can be described by an order parameter, 

M;i — TOiei cos(q- r) +7712^2 sin(q- r) +771363 + h.h., (1) 

where {ii} form an orthonormal triad and h.h. denotes 
"higher harmonics" such as terms proportional to sines and 
cosines of (2n + l)q • r with integer n. When the pitch vector, 
q, is normal to the plane of the rotating components, the rotat- 
ing components form a conventional helix''. ForTO3 = Osuch 
a state, which we call an 'ordinary helix' state, is observed in 
many rare-earth metals^", e.g. MnSi^' ~^, and FeGe^^'. We call 
a helix state with 7773 ^ 0, which is observed in some heavy 
rare-earth metals^*', a 'conical helix' state because the tip of 
the magnetization falls on the edge of a cone. A more compli- 
cated modulation arises when q lies in the plane of the rotating 
components. For 7713 = 0, we call such a state an 'ordinary 
cycloid' state because the profile of the magnetization resem- 
bles the shape of a cycloid. The state with 7713 7^ is called 
a 'conical cycloid' state. It is easy to see that the helical, but 
not the cycloidal, modulation preserves a residual symmetry 
under translations and suitable simultaneous rotations about 
the pitch vector 

Since M and P respectively break time reversal and spatial 
inversion symmetry, the leading P-dependent piece in a GL 
Hamiltonian density, hp, for a centrosymmetric, time reversal 
invariant system with cubic symmetry is^, 

ftp = P V2x + aP • M X V X M, (2) 

where x > a are coupling constants. We assume that 
P is a slave of M, in the sense that a non-zero P only oc- 
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curs due to the spontaneous development of a magnetic state 
with a non-zero M x V x M, which then, through the Hnear 
coupHng to P in (|2]i, induces a non-zero P. For an order pa- 
rameter ansatz given by Eq. [T] the macroscopic polarization, 
P, is given by minimizing the Hamiltonian density (|2]l over 
P, P = xamir7i2[e3 x q]. So P is normal to both q and the 
axis of rotation, 63. Note that in a conventional spin density 
wave state (mi or ni2 — 0), as in the helix states, P is zero. 
However, for a cycloid state, q _L 63, so there is a non-zero 
P. Note that P is entirely due to the cycloidal components 
nil and m2, and is independent of the uniform magnetization 
1113. Thus, while it is conceivable that magnetic fields strong 
enough to 'flop' the spins and the axis of rotation of the cy- 
cloidal components will alter P^'^'"^i^ , no explanation of how 
tuning the uniform component of M can affect the induced 
polarization has been offered. We will do so later in this pa- 
per. 

The paper is organized as follows: Section II lays out 
the Ginzburg-Landau Hamiltonian and the parameter regions 
which exhibits the cycloidal phase. Section III and V are de- 
voted to the phase diagrams of ordinary cycloidal state and 
conical cycloidal state respectively. In Section V, we explain 
why the reversal of the direction of the uniform magnetization 
in the conical cycloidal state can lead to the reversal of electric 
polarization. Section VI consists of conclusions. 



II. GINZBURG-LANDAU HAMILTONIAN 

We consider a Hamiltonian that is completely invariant un- 
der simultaneous rotations of positions and magnetization. 
This guarantees that any phase that can occur in our model 
is necessarily allowed in a crystal of any symmetry. The full 
Hamiltonian is given by, H = j {Hm + hp)dr = j hdv. Us- 
ing P = — xaM x V x M to eliminate P, we can write the 
total Hamiltonian density h entirely in terms of M, 

h = tm^ + uM^ + ifo (V • M)^ + ifi (V X M)^ 
+K2'M^ (V • M)^ + X3 (M • V X M)^ 
+Ki |M X V X M|^ 

+Dl\V {V ■M.)\^ + Dt\V {V xM)\'^, (3) 

where we have u, Dl,t > for stability. In Eq.[3] where the 
Landau expansion of the free energy is truncated at the fourth 
order, the usual gradient-squared term, c|VM|^, is omitted 
since, |VM|^ = (V • M)2 + |V X M|2, plus an unimpor- 
tant surface term which can be neglected. Notice that, for 
Kq = Ki and K2 = i^3 = K^, h is rotationally invariant 
in the spin space alone, so the Kis themselves are not pro- 
portional to the spin-orbit coupling constant (for e.g., via the 
above identity, i^Oi Ki ^ c). However, the difference among 
the A'j's should be small due to the smallness of the spin-orbit 
coupling. The effects of the competing magnetic interactions, 
which are present in the multiferroics and are responsible for 
the spatial modulation of M^-^ '^''^, are embodied in Kq, Ki, 
which can be negative leading to a spatially modulated order 
parameter For decoupled spin and coordinate spaces (Ki's 



equal), the energies of the helical and the cycloidal modula- 
tions of the spins are identical. In a system where the spin 
anisotropy constrains the spins to lie on a plane, and the lat- 
tice anisotropy forces q to be also on that plane, the energy of 
the cycloidal modulation can be lower than that of the helical 
modulation^. Such anisotropies have been implicitly taken 
as the driving force behind the cycloidal order by Mostovoy^, 
and Katsura et al}^. For cubic crystals, however, no such 
anisotropy exists among the principal directions. We argue 
below that, in this case, the magnetoelectric couplings them- 
selves, leading to the difference among the if/s, can lower 
the energy of the cycloidal state than that of any other state 
with an arbitrary angle between q and the plane of the mag- 
netization. 

Rather than exploring the complete parameter space of this 
model, we limit ourselves to two different parameter regions, 
which exhibit all the phases described above: 

Region I: Xo,i < 0, ^j>i small, t > , and 

Region II: t < 0, < 0, i^i > 0, K2 = = 0. 

We have checked that our results are robust against allow- 
ing small non-zero values of the various if^'s that we take to 
be zero. In that sense our results, in particular the topology 
of the phase diagrams shown in Figs. [T^ and|2^ for Regions I 
and II, respectively, and the orders of the various phase tran- 
sitions that we predict, are generic. As usual, our theoretical 
phase diagrams can be related to experimental ones by noting 
that all of the phenomenological parameters (i, Ki, Dl t, u) 
in our model should depend on experimental parameters like, 
e.g., temperature (T). Thus, an experiment in which, e.g., T is 
varied with all other parameters held fixed will map out a lo- 
cus of points through our theoretical phase diagrams. In Lan- 
dau theories, t is expected to vary from large positive values, 
corresponding to disordered phases with M(r) — 0, at high 
T, to smaller values at which M(r) 7^ become possible. In 
order to access the conical cycloid state, we must also allow 
Ko (T) and Ki (T) to change sign as T is decreased. 

For the most part we will work in mean field theory, which 
is simply finding a magnetization configuration M(r) that 
minimizes the Hamiltonian (O. Clearly, the task of finding 
the global minimum is a formidable one. Instead, we restrict 
ourselves to ansatzes of the form: 

M = TOiei cos(q • r) + 1112(22 sin(q • r) + Mq, (4) 

where the spatially constant vector Mq is allowed to point in 
any direction. (Given the global rotation invariance under si- 
multaneous rotations of magnetization and space, an infinity 
of other solutions trivially related to (|4]l by such rotations, and 
with exactly the same energy, also exist, of course.) In the spe- 
cial case of q along x direction (or, equivalently, anywhere in 
the X — y-plane), this is a cycloid state with a uniform back- 
ground magnetization Mq = {Mqi, Mq2, Mq^). When q is 
along z direction, it is a helix state. Inserting this ansatz (|4]l 
into the Hamiltonian ([3]), and integrating over the volume V of 
the system, we can obtain the energy of the system. Through 
the minimization of the energy, we find the conical cycloid 
state is the only state with a non-zero Mq when K3 < K4. 
In addition, the optimal direction for q is always either in the 
(x — y) plane, or orthogonal to it. Putting these facts together 
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means that all of the minimum energy configurations are of 
the form ([TJ. Furthermore, when q lies in the (x — y) plane, 
we can always use the global rotation invariance of our model 
to rotate q to lie along the a;-axis, and will henceforth do so. 



III. ORDINARY CYCLOID STATE 

In Region I, the dominant terms in the Hamiltonian involv- 
ing the uniform component are trri^ + um|, therefore the 
lowest energy states have = 0. Small negative Ki^i 
clearly cannot change this fact. The energy for the ordinary 
cycloid (OC) state is obtained by inserting ([TJ with = 
into the Hamiltonian 

E/V = Tl (q) ml + Tt (q) mj + u$(m?, m^), (5) 



where (g) = {t + hW + DLq^) /2, Tt [q] = 
{t + Kiq^ + Drq"^) /2, and $(m^,m^) 
3 (mf + m|) /8 + m\ra\/A. In writing this, we have 
neglected the higher harmonics in Eq. ([T]i, whose amplitude 
vanishes much faster (specifically, as fast or faster than |TOi|^) 
than the magnitude of the order parameter itself, and thus 
have negligible effects on the phase boundaries. For large 
positive t, all the terms in this energy are positive, and, 
hence, the lowest energy state is mi = m2 = 0; i.e., the 
paramagnet. As T decreases, t becomes smaller and the 
first phase transition that will occur depends on whether the 
minimum over g of {q) or Tt (q) becomes negative first. 
For r = Ki/Kq < ^/Dl/Dt, (g) becomes negative 
first at toLS — Kq/'^Dl, and nii starts to be nonzero. This 
boundary between paramagnet and the ordinary longitudinal 
spin density wave (OLS) phase (m2 = = 0, rni ^ 0) is 
the horizontal (solid blue ) line in the phase diagram Fig. [T^ 
in the (r, t) plane for fixed negative Ko and all Kiyi = 0. 

The OLS phase will, as we continue lowering t, eventually 
become unstable to a non-zero m2', this is the OC state. By 
minimizing the energy (|5]i in the OLS phase, we find q^ ~ 
ll.min = -Kq/2Dl and m\ = 2{toLS ~ t)/3u. Inserting 
these into (|5]l we find that the coefficient of m| becomes neg- 
ative below tLoc = toLS [3r - (1 + 3Dt/Dl) /2]. This 
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FIG. 1 : (Color online) (a) Phase diagram in Region I for the ordi- 
nary cycloid state. Solid (blue) lines represent second order phase 
transitions. Dotted (green) line indicates the first order transition 
to the helix state. The dotted (red) arrow represents one possible 
schematic locus of the experimental points obtained by varying T. 
To = (1 + 3Dt/Dl) /6. (b) The sequence of phases with decreas- 
ing T along the locus shown. 



value tioc of t therefore defines the locus of a continuous 
OLS-OC phase transition, and is the non-horizontal straight 
(solid (blue)) line in the r — t plane shown in Fig. [1^. 

For r > X Tt becomes non-zero first, which seems 

to imply that one enters the ordinary transverse spin density 
wave (OTS) phase (mi — ma = 0, m2 ^ 0) first for large 
r. However, it is not true because the OTS phase always has 
higher energy than the ordinary helical (OH) phase. The en- 
ergy for the ordinary helix state is 



E/V = TT{ml + ml) + $(m?, m^). 



(6) 



The minimization of the energy over the direction of (mi , m2 ) 
vector yields |mi| — \m2\ — m/f/V2, that is, a circular 
helix. Further minimization over and q gives the en- 
ergy EoH of the ordinary helix state Eqh/V = —{ton — 
t)2/4ufort < toH, where t OH = r'^DLtoLS / Dt- The 
energy for the OTS state is EqtsI^ = -{ton — t)'^/6u, 
which is obtained from equation ^ by setting mi = and 
1^ = iT.min = ~2st' ^^^^ minimizing over m2. Eqts 
is clearly higher than Eqh- Hence, the helical state is always 
favored over the OTS state throughout Region I of the phase 
diagram. Note that toH defines the boundary for the second 
order transition from the paramagnet to the OH state. 

There is also a direct first order phase transition between the 
OH and the OLS states along the line where Eqh = Eqls- 
Here Eols/V — ~{toLS — t)'^/Qu is the energy for OLS 
state obtained from equation (|5]l. This equality yields the first 

order phase boundary ioLH = (1/3/2*0^ — ioLs)/(\/3/2— 
1) between the OH and the OLS states (the dotted (green) 
line). The line for the OLS-OC transition always intersects 
the first order OLS-OH phase boundary before crossing the 
paramagnet-OLS boundary. This therefore always yields the 
topology shown in Fig. [T^. 

A typical experimental locus through this phase diagram, 
namely one in which t decreases as temperature T does, with 
r constant, is shown in Fig. [T^. The sequence of phases that 
results is illustrated in Fig. [TJ). We see that the paramagnet 
to ordinary cycloid phase transition is always preempted by 
a paramagnet to OLS phase transition, and the cycloid state 
is always elliptical. Both of these predictions are borne out 
by recent experiments on TbMnOg'^ '"'. On the other hand, a 
direct transition to the circular helix state is predicted by our 
theory, and has indeed been observed experimentally^i^^. 

All of the above statements are based on mean field the- 
ory, that is theory without considering the fluctuations. Go- 
ing beyond mean field theory, very general arguments due to 
Brazovskii-'* imply that, in rotation invariant models, any di- 
rect transition from a homogeneous state (paramagnet) to a 
translationally ordered one (OLS and OH) must be driven first 
order by fluctuations. Consideration of topological defects 
and orientational order^^'^^-^^ supports this conclusion, but 
raises the additional possibility that direct transition between 
the homogeneous and the translationally ordered phases could 
split into two, with an intermediate orientationally ordered 
phase, analogous to the 2D "hexatic" phase^*^. In the present 
context, this implies that both the paramagnet to OLS and OH 
phase transitions are either driven first order by fluctuations, or 
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split into two transitions with an intermediate orientationally 
ordered phase. Crystal symmetry breaking fields neglected in 
our model could invalidate this conclusion, if strong enough. 



IV. CONICAL CYCLOID STATE 
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In Region II, we can show that conical cycloid (CC) state 
of the form M = (toi cos{qx), m2 sm{qx), m^) is the low- 
est energy state among all the possible states with arbitrary 
mutual angles between the uniform magnetization, q, and the 
cycloid plane. The energy E for this state takes the form 



E/V = {t + Kaq^ + DLq^ + 2umi^mll2 
+ + Kiq^ + Drq 



(7) 



K:iq^ml) ml/2 



where we have again neglected the higher harmonics in Eq. 
([TJ. In this region, the h.h. terms do not vanish as the conical 
longitudinal spin density wave (CLS) (to2 = 0^ "^1,3, 7^ 0) or 
conical transverse spin density wave (CTS) (toi — 0, m2,3 ^ 
0) to FM transition in Fig. |2]is approached. However, we 
have verified that amplitudes of the h.h. terms are only a very 
small fraction of the cycloidal components mi and TO2 (not 
of the uniform component m^), therefore their neglect below 
(but close to) the lower cycloidal transition temperature of 26 
K is justified. They have little or no quantitative effect on our 
phase diagram or the orders of the transition. 

Since t < 0, we can minimize Eq. (Q over TO3 with 
TOi = 7712 = 0, and find a ferromagnetic (FM) state with 
ma = yj—t/2u. For large positive Kq and Ki, this fer- 
romagnetic state is clearly stable against the development of 
non-zero mi and m2. It also clearly becomes unstable against 
the development of a non-zero mi if K^) is lowered to nega- 
tive values, because then the coefficient {Koq^ + D^q^) of 
m\ becomes negative for sufficiently small q. This instability 
(which is clearly into the CLS state) will occur at = 0, at 
a wavevector q satisfying q\ ~ —Ko/2Dl- Note, how- 
ever, that now, because Kq is being varied through zero, this 
wavevector will now vanish as the transition is approached 
from below. The order parameter m^ = Kq/2uDl also 
vanishes as this transition is approached. Thus, this tran- 
sition is, like the (3 - incommensurate transition in quartz 
and berlinita^S, simultaneously a nucleation transition (q van- 
ishes), and an instability transition (order parameter vanishes). 
Indeed, this transition and the FM^ CTS transition, which 
is of the same type and will be discussed below, are, to our 
knowledge, the first examples of transitions that exhibit such 
a dual character in a model without terms linear in the gradient 
operator. 

We can find the loci of instability between the CLS phase 
and the CC state by calculating the coefficient of rri^ in 
O in the CLS phase, and finding where it becomes nega- 
tive. The minimization of the energy (|7]| over q, and mi 



yields q 



-Ko/2Dl, mi 



{t + Kl/2DL) /2u and 



m\ — Kq/2uDl. Inserting these expressions into (|7J and 
taking the coefficient of mj to be zero, we find the CLS to CC 



FIG. 2: (Color online) (a) Phase diagram in Region II for the conical 
cycloid state. Solid (blue) lines are the boundaries between different 
phases. The dotted (red) arrows represent possible paths for transi- 
tion to the CC state via continuous transitions, (b) The succession of 
the phases with decreasing T. The (green) arrow represents a direct 
first order transition between the FM and the CC state. 



phase boundary as: 



2 



Dt 



KqK, 
2uDl 



2u ' 



(8) 



Similar analysis of the sequence of the phase transition, FM 
CTS CC, yields the schematic phase diagram on the 
A'o — Ki plane given in Fig. |2^. The phase boundary be- 
tween FM and CTS is given by Ki = tKj,/2u. The phase 
boundary between the CTS and the CC phase at small Kq is 
Ki = 2Ko/ {Dl/Dt — 1), which is also shown in Fig.|2^. 

Fig. |2] shows that it is not possible to go from the FM to the 
CC state via a continuous transition, except at a single special 
point. Generic paths like the diagonal dashed lines in Fig. 
must go through either the CLS or the CTS state, so two tran- 
sitions are required to reach the CC state, which, additionally, 
must be elliptical. Hence the only way there can be a direct 
transition from the FM state to the CC state is via a first order 
phase transition, which is not addressed by our theory. This 
prediction is borne out by experiments of CoCr204, where the 
direct FM to CC transition is indeed first order—. 



V. MAGNETIC REVERSAL OF THE ELECTRIC 
POLARIZATION: 

The polarization P = xamim2y in the CC state is in the 
xy plane, normal to £3 and q. It is independent of the uniform 
magnetization, ma. Experimentally'^, the sample is cooled 
through Tc in the presence of a small electric field, E — Eoy, 
and a small magnetic field, H = Hqz. The direction of the 
pitch vector, x, or, equivalently, the axis of rotation, z, are 
set by the direction of P (E), which determines the 'helicity' 
of the cycloid^. It is found, at first, that P is uniquely deter- 
mined by E alone, independent of the initial direction of H, 
as expected. However, once P and m^ have set in, changing 
Ho to —Ho not only reverses the direction of ma, but also, 
quite unexpectedly, reverses the direction of P as well. In 
the literature^' this has lead to the definition of the 'toroidal 
moment', T = P x M, as the order parameter. 

It is clear that the experimental system is in the conical cy- 
cloid state, where ma, q and P are always in mutually orthog- 
onal directions'^. Further, as expected for this state, the di- 
rections of ma and P are uniquely determined by the small 
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FIG. 3: The reversal of the polarization (P) by the reversal of the 
magnetization (ma), (a) If ma rotates to — ms, remaining perpendic- 
ular to q, the cycloidal {xy) plane must rotate accordingly to always 
remain transverse to ms, which is the lowest energy configuration. 
Since P is in the cycloidal plane, it will rotate by a total angle tt. (b) 
An intermediate stage when ms has rotated by an angle |- and points 
in the y direction. At this stage, P points in the —z direction. 

cooling fields, H and E, respectively, which add terms to 
the Hamiltonian that split the degeneracy between the min- 
ima corresponding to the different directions. Now assume 
that the direction of H is reversed, iJp — > — ^fo> revers- 
ing the direction of 7713 once it has well developed. There 
are two ways the uniform magnetization can reverse its di- 
rection. First, m,3 may continue to remain along the z-axis 
and its magnitude may pass through zero to become —7773 for 
H — —Hqz. If this is the case, P will remain fixed in the di- 
rection y, since the mutual orthogonality of 7773, q and P can 
always be maintained and there is no direct coupling between 
7773 and P. However, since 7773 is already well developed and 
large (T,„ = 93 K), due to the magnetic exchange energy cost 
it may be energetically more favorable to leave the magni- 
tude of 777.3 unchanged, and its direction may rotate in space 
to —z. If this is the case, then 7773 must rotate staying on the 
y ~ z plane, since that way it always remains perpendicular to 



q, whose direction fluctuations cost the crystalline anisotropy 
energy. It is then clear, see Fig. [3] that the cycloid plane itself, 
which is always perpendicular to 7773 to maintain the lowest 
energy configuration, must rotate about x by a total angle tt. 
It follows that P, always on the cycloid plane, reverses its 
direction to —y. This way, even though there is no dynami- 
cal coupling between 7773 and P, the latter can also rotate by 
an angle tt as a result of the former reversing its direction in 
space. Based on this, we predict that, at some intermediate 
H ~ —H'z, where H' < Hq, P points in the direction — z, 
which can be experimentally tested. 



VI. CONCLUSIONS 

To conclude, we've shown that the magnetic cycloidal or- 
ders, and the resulting multiferroicity, can naturally arise due 
to the magnetoelectric couplings even in rotationally invari- 
ant systems, or in cubic crystals. This explains such orders in 
CoCr204, which lack easy plane anisotropics, and are hence 
outside the realm of the previous theoretical studies on multi- 
ferroics. We also predict that a second order transition from 
the ferromagnet to the conical cycloid state can only occur 
through an intervening conical longitudinal or transverse spin 
density wave state with the ultimate cycloidal state being el- 
liptical. A direct such transition, then, must be first order An 
important feature of our Ginzburg-Landau theory is that we 
do not need to invoke an arbitrary (and ad hoc) 'toroidal mo- 
ment' to explain the interplay between the magnetization and 
the polarization - the behavior which has been attributed to 
the toroidal moment arises naturally in our theory. 

We thank D. Drew, D. Belitz, and R. Valdes Aguilar for use- 
ful discussions. This work is supported by NSF, NRJ, LPS- 
NSA, and SWAN. 
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